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Abstract 

In this paper we present a model based on dynamics of the electrons in 
the plasma using a simplified Boltzmann equation coupled with a Poisson 
equation. 

The motivation arose to simulate active plasma resonance spectroscopy 
which is used for plasma diagnostic techniques, see [2], [14] and [16] . 

We are interested on designing splitting methods to the model prob- 
lem. 

First we reduce to a simplified transport equation and start to analyze 
the abstract Cauchy problem based on semi-groups. 

Second we extent to the coupled transport and kinetic model and apply 
the splitting ideas. 

The results are discussed with first numerical experiments to give dis- 
cuss the numerical methods. 

Keywords: kinetic model, neutron transport, dynamics of electrons, trans- 
port equation, splitting schemes, semi-group. 
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1 Introduction 

We motivate our studying on simulating a active plasma resonance spectroscopy, 
which is well established in plasma diagnostic techniques. 

To study the model with simulation models, we concentrate on an abstract 
kinetic model, which described the dynamics of electrons in the plasma by using 
a Boltzmann equation. The Boltzmann equation is coupled with the electric field 
and we obtain coupled partial differential equations. 
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Our combined model is done simplified to apply with functional analytical 
tools. We discuss the description of a positive semi-group, which helps to do 
the numerical estimations in the splitting schemes. 

Second a numerical method is discussed with respect to separate differential 
and integral part of the equations. 

The numerical approximation is done by applying splitting methods of sec- 
ond order. 

The paper is outlined as follows. 

In section [2] we present our mathematical model and a possible reduced 
model for the further approximations. 

The functional analytical setting with semi-groups are discussed in section[3J 
The splitting schemes are presented in in Section|4]and the numerical integration 
of the integro-part is discussed in Section [5] 

Numerical experiments are done in Section[6] In the contents, that are given 
in Section [3 we summarize our results. 



2 Mathematical Model 

In the following a model is presented due to the motivation in [2] , [H] and [16] . 

The models consider a fluid dynamical approach of the natural ability of 
plasmas to resonate in the near of the electron plasma frequency u pe . 

Here we specialize to an abstract kinetic model to describe the dynamics of 
the electrons in the plasma, that allows to do the resonance-analysis. 

The Boltzmann equation for the electron particles are given as 

= -v ■ V x f(x, v, t) - — • V„/(s, v, t) 
at m e 



-a(x,v,t)f(x,v,t)+ n(x,v,v')f(x, v',t) dv' , (1) 
Jv 

f(x,v,0) = f (x,v), (2) 



and boundary conditions are postulated at the boundaries of P (plasma). 

In front of the materials we assume complete reflection of the electrons due 
to the sheath f(vu + v±) with vu is the parallel and v± perpendicular to the 
surface normal vector. tf> is the electric field. 

The Boltzmann's equation has to be coupled with the electric field. The 
electrostatic approximation of the field is represented by the potential that is 
valid on the complete volume S. 

We apply the the Poisson's equation: 

-V a -(eW) = { <»i-ff*S) (3) 

the permittivity is equal to eo in the plasma P and eoe_D in the dielectric D. 
<fi fulfills the boundary conditions U n at any electrode E n and — n • V0 = at 
isolator /, whereas n is the normal vector of the isolator surface. 
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On the surface of the dielectric a surface charge a may accumulate and leads 
to a transition condition: 

A(eV<£) = -a. (4) 



3 Semi-groups for Transport Equations 

In the following, we derive the exponential growth of the transport semi-groups 
that is used in the section of the numerical methods. 

We discuss in the following subsections two directions of the transport regimes: 

• Neutron transport and 

• Electron transport. 



3.1 Transport model for the neutrons 

For this model we can assume that /(x, v, t) describe the density distribution of 
particles at position ieS with speed v £ V at time t G [0, T], see also [I] and 

El- 

The space S is assumed to be a compact and convex subset of IR 3 with 
nonempty interior, and the velocity space V is: 
V := {v e 1R 3 : v min < |||| 2 < v max } 
for v min > and v max < oo. 

Assumption 3.1. We have the following assumptions: 

• Particles move according to their speed v. 

• Particles are absorbed with function a (e.g. probability function), depend- 
ing on x and v. 

• Particles are scattered to a scattering kernel k depending on position x, 
incoming speed v' and outgoing speed v. 

The neutron transport is given as: 

JL ' = - v . V/(x, v, t) - a(x, v, t)f(x, v, t) 

+ / K(x,v,v')f(x,v',t)dv' 7 (5) 
Jv 

f(x,v,0) = f a (x,v), (6) 

and boundary conditions are included in the transport operator Aq see in the 
following abstract Cauchy problem. 

In the following we deal with the abstract Cauchy Problem for the simplified 
model. 
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3.1.1 Abstract Cauchy problem: Transport model for the neutrons 

We have a Banach space X := L 1 (S x V) with Lebesgue measure onSxl^C IR 6 
and define the abstract Cauchy problem as: 

= Bu(t) , (7) 

= (A - M a + K K )u(t) , (8) 
u , (9) 

where u £ X. 

We have the following operators: 

1. ) Collision- less transport operator 

2. ) Absorption operator 

3. ) Scattering Operator 

An important results for further numerical analysis is the fact, that the 
transport semi-group can be estimated by an exponential growth, see [3]: 

Corollary 3.1. We assume that s(B) > —oo is a dominant eigenvalue and 
(S(t))t>ois irreducible Then the transport semi-group (S(t)) t >o has balanced ex- 
ponential growth. There exists a one- dimensional projection P satisfying < Pf 
whenever < / such that: 

\\exp(-s(B)t)S(t) - P\\ <Mexp(-ei), (10) 
for all t > and appropriate M > 1 and e > 0. 

3.2 Transport model for the electrons or ions 

For this model we can assume that f(x, v, t) describe the density distribution of 
particles at position x S S with speed v G V at time t 6 [0,T], see also [4] and 

The space S is assumed to be a compact and convex subset of IR 3 with 
nonempty interior, and the velocity space V is: 
V := {v <E H 3 : v min < |||| 2 < v max } 
for v min > and v max < oo. 

Assumption 3.2. We have the following assumptions: 

• Particles move according to their speed v. 

• Particles are absorbed with function a (e.g. probability function), depend- 
ing on x and v. 

• Particles are scattered to a scattering kernel k depending on position x, 
incoming speed v' and outgoing speed v. 

• Particles are influenced by the static electric field 4>, which can be derived 
by the kinetic theory. 



dult) 

dt 
du(t) 

dt 
u(0) = 
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The electron transport is given as: 

df(x ' v ' t) = -v ■ V x f{x, v, t) - — v x <f> ■ V v f(x, v, t) 
ot m e 

—a(x,v,t)f(x,v,t)+ / K(x,v,v')f(x,v',t)dv\ (11) 
Jv 

f(x,v,0)=f (x,v), (12) 

and boundary conditions are included in the transport operators. (f> is the 
electric field. 

Further we have the Poisson's equation: 



-V,. W ) = { e ^-^ ^ 



(13) 



in£> ' 

the permittivity is equal to eo in the plasma P and e €o in the dielectric D. 

df{x ^' t] = -v ■ V x /(x, v, t) - V X D ■ V x f{x, t) 
ot 

—a(x,v,t)f(x,v,t)+ / K(x,v,v')f(x,v',t)dv', (14) 
Jv 

f(x,v,0) = f (x,v), (15) 

and boundary conditions are included in the transport operators and A\ see 
in the following abstract Cauchy problem. D is the diffusion parameter that 
includes the electric field. 

Next we deal with the abstract Cauchy Problem for the simplified model. 

3.2.1 Abstract Cauchy problem: Transport model for the neutrons 

We have a Banach space X := L 1 (S x V) with Lebesgue measure on S x V C H 6 
and define the abstract Cauchy problem as: 

du(t) 



dt 
du{t) 



= Bu{t) , (16) 
= (A +A 1 -M a + K K )u(t) , (17) 



dt 

u(0) - u , (18) 

where u G X. 

We have the following operators: 

1. ) Collisiondess transport operator 

2. ) Diffusion operator 

3. ) Absorption operator 

4. ) Scattering Operator 

An important results for further numerical analysis is the fact, that the 
transport semi-group can be estimated by an exponential growth. 
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Corollary 3.2. We assume that s(B) > — oo is a dominant eigenvalue and 
(S(t)) t >ois irreducible Then the transport semi-group (S(t)) t >a has balanced ex- 
ponential growth. There exists a one- dimensional projection P satisfying < Pf 
whenever < / such that: 

\\exp(-s(B)t)S(t) - P\\ <Mexp(-et), (19) 

for all t > and appropriate M > 1 and e > 0. 

In the next section we discuss the splitting schemes. 



4 Splitting schemes 

The operator-splitting methods are used to solve complex models in the geo- 
physical and environmental physics, they are developed and applied in [T5], [20J 
and [21 j . This ideas based in this article are solving simpler equations with 
respect to receive higher order discretization methods for the remain equations. 
For this aim we use the operator-splitting method and decouple the equation as 
follows described. 



4.1 Splitting methods of first order for linear equations 

First we describe the simplest operator-splitting, which is called fi for the fol- 
lowing system of ordinary linear differential equations: 

d t c(t) = A c(t) + B c(t) , (20) 

whereby the initial-conditions are c" = c(t n ). The operators A and B are 
spatially discretized operators, e.g. they correspond to the discretized in space 
convection and diffusion operators (matrices). Hence, they can be considered 
as bounded operators. 

The sequential operator-splitting method is introduced as a method which 
solve the two sub-problems sequentially, where the different sub-problems are 
connected via the initial conditions. This means that we replace the original 
problem (f2"U|) with the sub-problems 

dc* (t) 

' [ '=Ac*{t), with c*(t n ) = c n , (21) 



dt 

8r**(t\ 

[ > = Bc**(t) , withc**(r) = c*{t n+1 ) , 



dt 

whereby the splitting time-step is defined as r„ = t n+1 — t n . The approximated 
split solution is defined as c n+1 = c**(t n+1 ). 

Clearly, the change of the original problems with the sub-problems usually 
results some error, called splitting error. Obviously, the splitting error of the fi 
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method can be derived as follows (cf. e.g.|I2]) 

Pn = -(exp(T n (A + B)) - exp(r n B) exp(r„A)) c(t n ) 

T 

= l -T n [A,B]c{t n ) + 0{T 2 ) . (22) 

whereby [A, B] := AB — BA is the commutator of A and B. Consequently, the 
splitting error is 0(r n ) when the operators A and B do not commute, otherwise 
the method is exact. Hence, by definition, the fi is called first order splitting 
method . 



4.2 Sequential splitting method for non-linear problems 

We could use the result for the general formulation of nonlinear ordinary differ- 
ential equations: 

c'(t) = F 1 (t,c(t)) + F 2 (t,c(t)) , (23) 

where the initial-conditions are given as c™ = c(t"). 

As before, we can decouple the above problem into two (typically simpler) 
sub-problems, namely 

Be* (t) 

_ii = Ft(t, c*{t)) with t n < t < t n+1 and c*{t n ) = c n , (24) 
F 2 (t,c**(t)) withi™ < t < t n+1 andc**(t") = c*(t n+1 ) ,(25) 



dt 

where the initial- values are given as c" = c{t n ) and the split approximation on 
the next time level is defined as c n+1 = c**(t n+1 ). 

For this case the splitting error can be defined by use of the Jacobians of the 
non- linear mappings F\ and F 2 , namely as 

1 BE BE 

Pn = z^F,, -^- Fl ]{t n ,e{n) + 0(jl) . (26) 
Hence, for the general case the splitting error has of first order, i.e. 0(r„). 

4.3 Higher order splitting methods for linear operators 

So far we defined the fi which has first order accuracy. However in the practical 
computations in many cases we require splittings of higher order accuracy. 

4.3.1 Symmetrically weighted sequential splitting. 

In the following we introduce a weighted sequential splitting method, which is 
based on two sequential splitting methods with different ordering of the op- 
erators. I.e. we consider again the Cauchy problem (|20[) and we define the 
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operator-splitting on the time interval [t n ,t n+1 ] (where t n+1 = t n + t„) as fol- 
lows 



dc*(t) 

dt 
dc**{t) 
dt 



and 



dv*{t) 

dt 
dv**(t) 

at 



Ac*(t) , withc*(f l ) =c n , (27) 
= Bc**(i) , with c**(t n ) = c*{t n+l ) . 

= Bv*(t) , with v*(t n ) = c" , (28) 
= Av**(t) , with v**{t n ) = v*{t n+1 ) . 



where c™ is known. 

Then the approximation at the next time-level t n+1 is defined as 

„** (j-n+l\ I ** (-in+l\ 
c n+l = c_Jt )+ V V L (29) 

The splitting error of this operator splitting method is derived as follows (cf. 

®) 

1 1 
Pn = — {exp(r„(A + B)) - - [exp(r n B) exp(r„^) + exp(r„A) exp(r„B)]} c(t n ) 

Tn 2 

= 0(r 2 ). (30) 

An easy computation shows that in general case the splitting error of this 
method is 0(r 2 ), i.e. the method is of second order accurate. (We note that in 
the case of commuting operators A and B the method is exact, i.e. the splitting 
error vanishes.) 

4.3.2 Strang-Marchuk splitting method. 

One of the most popular and widely used operator-splittings is the so-called 
Strang splitting (or Strang-Marchuk splitting), defined as follows |15i 119] - The 
methods reads as follows 

W = Ac*{t) , with t n <t< i n+1 / 2 and c*(t n ) = c™ , (31) 



dt 

r)r** (t) 

W = Bc **( t ) t with t n <t< t n+1 and c**(t n ) = c*(t n+1/2 ) 



dt 

gt K ' = Ac***(t) , with t n+1 ' 2 < t < t n+1 and c ***{t n+1 ' 2 ) = c**(t n+1 ) , 

where i' i+1 / 2 = t n + 0.5t„ and the approximation on the next time level t n+1 is 
defined as c' i+1 = c***{t n+1 ). 
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The splitting error of the Strang splitting is 

Pn = ^(r n ) 2 ([B, [B, A]} 2[A, [A, B]]) c(t n ) + 0(r„ 4 ) . (32) 

(See, e.g. ( [10 . ) This means that this operator-splitting is of second order, 
too. (We note that under some special conditions for the operators A and B, 
the Strang splitting has third order accuracy and even can be exact [5]. ) 

In our application the first order splitting for the convection-reaction- and 
the diffusion-dispersion-term are applied, because of the dominance of the space- 
error. The time-error for this combination was only a constant in the total error. 

In the next subsection we present the iterative-splitting method. 



4.4 Iterative splitting method 

The following algorithm is based on the iteration with fixed splitting discretiza- 
tion step-size r, namely, on the time interval [t n ,t n+1 ] we solve the following 
sub-problems consecutively for i = 0,2, . . . 2m. (Cf. [T3] and [5].) 



^> = Aa(t) + Bci_i(t), with Ci (t n ) = c n (33) 
and c (t n ) = c" , c_! = 0.0, 

?*±M= Aci (t) + B* +1 (t), (34) 
with Ci +1 {t n ) =c" , 

where c" is the known split approximation at the time level t — t n . The split 
approximation at the time-level t — t n+1 is defined as c n+1 = c 2m +i(t n+1 ). 
(Clearly, the function Ci + i(t) depends on the interval [t n ,t n+1 ], too, but, for the 
sake of simplicity, in our notation we omit the dependence on n.) 

In the following we will analyze the convergence and the rate of the conver- 
gence of the method (|3"3")) - l|3"4"|) for m tends to infinity for the linear operators 
A, B :X — > X where we assume that these operators and their sum are genera- 
tors of the Co semi-groups. We emphasize that these operators aren't necessarily 
bounded, so, the convergence is examined in general Banach space setting. 

Theorem 4.1. Let us consider the abstract Cauchy problem in a Banach space 
X 

d t c{t) = Ac[t) + Bc{t), 0<t<T 
c(0) = c o 

where A, B, A+B :X — > X are given linear operators being generators of the Co- 
semi-group and Co G X is a given element. Then the iteration process iS3\) - (3^ 
is convergent and the and the rate of the convergence is of second order. 
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Proof. Let us consider the iteration (f3"3")) - ([3~4l on the sub-interval [t n ,t n+1 ]. For 
the error function ej (t) = c(t) — Ci (t) we have the relations 



e t (t n )=0 



and 



(36) 
(37) 



d t e i+1 (t) = Aei(t) + Be i+1 (t), te(t n ,t n+1 ], 
e i+1 (t n ) = 

for m = 0,2,4,..., with eo(0) = and e_i(t) = c(t). In the following we 
use the notations X 2 for the product space X x X enabled with the norm 
||(«,u)|| = max{||u||, ||u||} (u,v e X). The elements £i{t), Ji(t) E X 2 and the 
linear operator A : X 2 — > X 2 are defined as follows 



£i(t) = 



ei(t) 
e»+i(*) 



; ^i(t) = 



ei-i(t) 




4 = 



A 
A B 



(38) 



Then, using the notations (f3"5|). the relations (|78p - (j3"7| can be written in the 
form 

d t £ i {t) = A£ i {t)+F i {t), t£(t n ,t n+1 ], 
£i(t n ) = 0. 

Due to our assumptions, A is a generator of the one-parameter Co semi-group 
(A(t)) t >o, hence using the variations of constants formula, the solution of the 
abstract Cauchy problem (|39[) with homogeneous initial condition can be written 

as t 

Si(t) = j exp(.4(i - a))^i{a)da, t e [t n ,t n+1 ], (40) 

(See, e.g. [4j.) Hence, using the denotation 

Halloo = sup te[t „ !t „ +1] ||^(t)|| (41) 

we have ^ 

\\£i\\(t)< \\F i \\ 00 Jjexp(A(t-s))\\ds = 

(42) 



= ]|e i _ 1 || yj|exp(^(t- S ))l|d S , t€[t n ,t n+1 }. 
Since (A(t)) t > is a semi-group therefore the so called growth estimation 

|| exp(.At)|| < Kexp(ut); t > 
holds with some numbers K > and cj e IR [4J. 



(43) 



• Assume that (A(t)) t >o is a bounded or exponentially stable semi-group, 
i.e. (|43|) holds with some ui < 0. Then obviously the estimate 



exp(„4i)|| < K; t>0 



(44) 
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holds, and, hence on base of ([75)1 , we have the relation 

\\£i\\(t) < Krje^l te(0,T„). (45) 



• Assume that (A(t)) t >o has an exponential growth with some ui > 0. Using 
(|78]) we have 

||exp(>l(t-«))||d« te[f,t" +1 ], (46) 



^(t) = -(exp(u;(i-f"))-l), te[f,f +1 ]. (47) 
< - (exp(wr„) - 1) = Kr n + 0(t*) (48) 



where 



Hence 



The estimations (|45|) and (1481) result in that 

\\£i\\oo=KT n \\e i - 1 \\ + 0(T*). (49) 
Taking into the account the definition of £i and the norm |j • Hoc, we obtain 

\\e i \\=Kr n \\e^ 1 \\+0{Tl), (50) 

and hence 

||e i+1 || =K 1 rl\\e l ^\\+0{rl), (51) 
which proves our statement. 

□ 

Remark 4.1. When A and B are matrices (i.e. (|33[) -(|34 |) is a system of the 
ordinary differential equations), for the growth estimation (|43|) we can use the 
concept of the logarithmic norm. (See e.g.|llj.) Hence, for many important 
class of matrices we can prove the validity of (l43l with lo < 0. 



Remark 4.2. We note that a huge class of important differential operators 
generate contractive semi-group. This means that for such problems -assuming 
the exact solvability of the split sub-problems- the iterative splitting method is 
convergent in second order to the exact solution. 

Modify to : 
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d ° 1 ^ - Ac\t) + Bc 1 - 1 ^), with c*(t n ) = c i - 1 (t n+1 ) (52) 
rting values c°(t n ) = c(t n ) results of last iteration , c~ 1 (t") = 0.0, 
Ac\t) + Bc l+1 {t), (53) 



dt 

and the starting values c (t n ) = c(t") results of last iteration , c~ 1 (t") = 0.0 

d c ' +1 (*) = 4r*m 4- R^ +1 
Si 

with c i+1 (t n ) = c i (t n+1 ) , 

e > \c l+1 {t n+1 ) - c^^r+^IStop criterion (54) 
result for the next time-step (55) 
c(t n+1 ) = c m (t n+1 ), for m fulfill the stop-criterion (56) 

for each i — 0,2, ... , where c™ is the known split approximation at the previous 
time level. 

5 Numerical Integration of the Integro-Part 

We deal with the following intcgro-differcntial equation: 



/ u(s) ds, (57) 
Jo 



du 
dt 

u{0)=u , (58) 



The integration part is done numerically with: 
Trapezoidal rule: 



/' 

J a 



k=l 



where the sub-intervals have the form [kh,(k + l)h], with h = (ba)/n and 
k = 0,l,2,...,nl. 

The higher order formulas are given as closed NewtonCotes formulas are 
given as 

where /j is a shorthand for f(xi), with Xi = a + i(b~a)/n, and n the degree. 
We obtain the following formulas for the Trapezoidal-rule: 

— =t/2(u(0) + u(t))ds, (60) 
«(0) - uo, (61) 

and obtain the analytical result: 

u(t) = |ex P (^) U (0)-^«(0), (62) 
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Degree 
1 
2 
3 
4 


Common name 
Trapezoid rule 
Simpson's rule 
Simpson's 3/8 rule 
Boole's rule 


Formula 

V(/o+/i) 
^(/o + 4/i + / 2 ) 

V(/0 + 3/!+3/ 2 + /3) 

^f(7/ + 32/ 1 + 12/ 2 
+32/ 3 + 7/ 4 ) 


Error term 

( b -f /(a)(0 

(6-a) 5 ,(4)/- 

2880 j w 

(b-a) 5 f(4)/^ 
6480 ■> v«J 
(b-a) f(6)(c\ 
1935360 J VS.; 



Table 1: Numerical Integration formulas (Closed NewtonCotes Formulas). 



For the higher order formula like Simpsons-rule, we have the following re- 
sults: 



du 
~dt 

U(0) = U 



t/Q(u(0) + Au{t/2) + u(t)) ds, 



(63) 
(64) 



We apply the idea of the polynomial solution: 
u(t) — a Q + a x t + a 2 t 2 + a 3 t 3 + . .. 

and we obtain the results with deriving the coefficients : 

di + 2a 2 t + 3a 3 t 2 + . . . 
= t/6(a + 4(a + + a 2 t 2 / 4 + a 3 t 3 / 

+a + a x t + a 2 t 2 + a 3 t 3 + . . . ) , 

a = Uq, 

and we obtain via coefficient comparison: 



■ + ...) 



= as = an = . . . = 



a = U 
ai = a 3 = a 5 
a 2 = 3a 
1 

«4 = j^a 2 , 



(65) 



(66) 



(67) 
(68) 
(69) 

(70) 



Remark 5.1. Such fast algorithms of generalized Taylor series about a function 
(here we apply numerical integration formulas) are computed very efficient, see 
also the decomposition ideas of JI]/. 



6 Experiments for the Plasma resonance spec- 
troscopy 

6.1 First Example: Matrix problem with integral term 

We deal with a simpler integro-differential equations: 
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c'{t) = c+ / c{t)dt, t g [0,1], (71) 



where we assume J c(t)dt = tc(t) as a first order approximation of the integral 
and deal with: 

c'(t) = c + tc(t), t G [0,T], (72) 
c(0) = c(to) = 1, (73) 

where T = 10.0 and we have the analytical solution for the approximation which 
is given as: 

t 2 

c(t) = exp(t+-X0), (74) 

We split into: 

A = 1, (75) 
B(t) = f, (76) 

We have the following solutions for the iterative scheme: 

Cl (t) = exp(A(t n+1 - t))c(t n ), i G (t n ,t n+1 ], (77) 



ca(t) =exp(/; n B(t)dt)c{t n ) 

+ J t !T +1 exp(/ s 4 " +1 S(i" +1 - t))tftAci(*)ds, * G (t n ,t n+1 ]. 



(78) 



it" 



where n — 0,1,..., N and t — T while the time-steps are given as At = 

f n+l _ t n_ 

We deal with the following recurrence relations with even and odd iterations: 
for the odd iterations: i = 2m + 1, 
for to = 0, 1, 2, . . . 

d(t) =exp(A(t-n)c(n 

+ f* n cxp(sA)B(s)c^ 1 {t n+1 -s)ds, te(t n ,t n+1 }. [ ' 

For the even iterations: i = 2to, 
for to = 1,2,... 

d(t) =exp(J tn B{s)ds)c(t n ) , gQ , 
+ J t ieMlo B{t)dt)Ac^ 1 (t n+1 -s) ds, te{t n ,t n+1 }. 

In the table [2] we obtain the numerical results of the iterative splitting 
scheme. 

In the Figure [U we present the one-side and two-side iterative results. 

Remark 6.1. In the experiments, we obtain improved results with each addi- 
tional step. By the way the solution blows up and we have to use also very fine 
time-steps to control the errors. Optimal results are obtain by using the integral 
part (stiff part) as the implicit part in the iteration (one-side over B). 
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At=l 


At=0.5 


At=0.25 


At=2- 3 


At=2~ 4 




1.7634 


0.4958 


0.1793 


0.0753 


0.0343 


C2 


0.8628 


0.1444 


0.0282 


0.0061 


0.0014 


C3 


0.2220 


0.0104 


5.2127e-04 


2.8455e-05 


1.6511e-06 


c 4 


0.1116 


0.0041 


1.8660e-04 


9.7846e-06 


5.5769e-07 


C5 


0.0971 


0.0039 


1.8367e-04 


9.7418c-06 


5.5644c-07 


C6 


0.0956 


0.0039 


1.8365c-04 


9.7418e-06 


5.5644c-07 


C7 


0.0954 


0.0039 


1.8365c-04 


9.7416c-06 


5.5644c-07 


C8 


0.0954 


0.0039 


1.8366e-04 


9.7343e-06 


5.5556e-07 



Table 2: Numerical experiment with 10 iterative steps for the first example. 



6.2 Real-life problem 

In the following subsections, we present our experiments based on the neutron 
transport. 

A simplified one-dimensional model is given as: 

8 t c + v8 x c-D8 xx c + ac=[ K{ x,v,v')c(x, V ', t )dv', 



The velocity v and the diffusion D is given by the plasma model. The initial 
conditions are given by c(x, 0) = c (x) and the boundary conditions are trivial 
d n c(x,t) = 0. 

A first integral operator is given as: 

/ k(x, v, v')c(x, v',t) dv' = / c(x,t)dt, 
Jn Jo 



A second integral operator is given as: 

We assume a simple collision operator: k(x, v, v') = q{v')(l + v' 2 ) 
where q(v') is the potential, e.g.v' 2 . 

We deal with the first integral operator and define the following operators: 
A = vU l - l V-D^[l-2l]I 
B = {-a + t)I 
while 

cxp(Bt) = cxp((-crt + t 2 /2)I) 

where / is the identity matrix of rank. 



6.2.1 One phase example 

The next example is a simplified real-life problem for a neutron transport equa- 
tion, which includes the gain and loss of a neutron 
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We concentrate on the computational benefits of a fast computation of the 
iterative scheme, given with matrix exponentials. 

The equation is given as: 



d t c + V-Fc 



-Aic 



\2c(x,t)dt, in CI x [0,t], 



F = v DV, 

c(x,£) = co(x), on Cl, 

c(x, i) = ci (x, i) , on dCl x [0, i] , 



(81) 

(82) 
(83) 
(84) 



In the following we deal with the semi-discretized equation given with the 
matrices: 



d t C = (A-Ai+A 2 )C, 



(85) 



where C = (ci,... , c/) T is the solution of the species in the mobile phase in 
each spatial discretization point (i = 1, . . . , I). 

We have the following two operators for the splitting method: 



D 



/ -2 1 
1 -2 



V 

/ 1 

-1 1 



1 -2 1 
1 -2/ 



+ 



v 

Ax 



\ 



V 



-1 1 

-1 1 J 



e IR 



7x7 



where I is the number of spatial points. 

/ Ai 
Ai 



Ai = 



V 



Ai 
Xi J 



7x7 



For the integral term we have the following ideas: 
Case 1: 

Jq \2c(x,t)dt w \2tc(x, t) 



(86) 



(87) 



(88) 
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and we obtain the Matrix: 

( X 2 t 2 /2 

\ 2 t 2 /2 



Aa 



V 



e JR 



Ixl 



(89) 



X 2 t 2 /2 

A 2 i 2 /2 J 

For the operator splitting scheme we apply A and B = — Ai + A 2 and we 
apply the iterative splitting method, given in equations ([75])- ([SO]). 
Case 2: 

We integrate the operator B with respect to the previous solutions Cj_i and 
we obtain the Matrix: 

A 2 (C,_ 1 ) 

J * A 2 ci ij; _i(x, s) ds ... ' 

/ A 2 c 2 ,i-i(x, s) ds 



V 







Jo A 2 c/,j_i(a; ) s) y 
(90) 



We obtain S(C) = A 2 (C < _i) + AiC 
The iterative scheme is given as: 
For i = 1,2,... 

d(t) = exp(A{t ~ t n ))C{t n ) 



+ //„ exp((i - a)A)B(Ci-i(a)) ds, t e [t n ,t n+1 ] 



(91) 



For the reference solution, we apply a fine time- and spatial scale without 
decoupling the equations. 

The Figure [5] present the numerical errors between the exact and the numer- 
ical solution. Here we obtain optimal results for one-side iterative schemes on 
operator B, means we iterate with respect to B and use A as right hand side. 

Remark 6.2. For all iterative schemes, we can reach faster results as for the 
The iterative schemes with fast computations of the exponential matrices stan- 
dard schemes. With 4 — 5 iterative steps we obtain more accurate results as 
we did for the expensive standard schemes. With one-side iterative schemes we 
reach the best convergence results. 



7 Conclusions and Discussions 

We present the coupled model for a transport model for deposition species in a 
plasma environment. We assume the flow field is computed by the plasma model 
and the transport of the deposition species with a transport-reaction model. 

Such a first model can help to understand the important modeling of the 
plasma environment in a CVD reactor. 
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Figure 1: Numerical errors of the one-side and two-side Splitting scheme: one- 
side splitting over A (upper figure), one-side splitting over B (middle figure) 
and two-side splitting scheme alternating between A and B (lower figure) with 
f , . . . , 8 iterative steps. 
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Figure 2: Numerical errors of the one-side Splitting scheme with A (upper 
figure), the one-side Splitting scheme with B (middle figure) and the iterative 
schemes with 1, . . . , 6 iterative steps (lower figure). 



